The Chemical Evolution of Galaxies 


[Tinsley 1980, Fund. Cosmic Phys., 5, 287] 
[Rana 1991, ARAA, 29, 129] 


There a number of processes that affect the chemistry of galaxies. 
These include 


e Type II Supernovae: In addition to creating iron-peak, and 
r-process elements, supernovae from high-mass stars produce the 
majority of a galaxy’s a-process elements (oxygen, sulfur, neon, 
argon, etc.) 


e Type Ia Supernovae: In the collapse of a white dwarf, virtu- 
ally the entire star will be changed into iron-peak elements. It is 
generally assumed that most of the iron seen in a galaxy originated 
in these objects. 


e Novae: Nova outbursts eject the products of CNO burning 
directly into space, hence they are a major contributer of nitrogen. 
Novae with high-mass white dwarfs can also produce a significant 
amount of neon, aluminum, magnesium, and sodium. 





e AGB stars: Some mixing does occur in the thermal-pulsing 
phase of AGB stars. This will enhance the CNO abundances of a 
system (especially nitrogen) and can build-up s-process elements. 





Because each these processes involves a different timescale, the 
elemental ratios (such as [O/Fe]) will be population dependent. 
Models for the build-up of these elements are more complex that 
those for a galaxy’s photometric evolution. However, they can 
still be written in closed forms via a series of coupled differential 
equa- tions. Furthermore, with a few approximations, analytical 
solutions to chemical evolution can be found. 


Metallicity Indicators 


There are only a few ways to measure (or estimate) the metallicity 


of a a galaxy: 


e Spectroscopy of Bright Stars: 


Advantages: 


Disadvantages: 


Measures chemical history of ISM via stars 
Good for both Fe-peak and a-process elements 
Requires high S/N absorption line spectra 
Only possible in very nearest galaxies 

Only possible for brightest (non-LTE) O-stars 


e Spectroscopy of H II Regions: 


Advantages: 


Disadvantages: 


Requires moderate S/N of emission-lines 
Sources are easy to find and measure 
Analysis is straightforward 

Only measures present day ISM 

Only possible in star-forming systems 
Mostly limited to a-elements (best for O) 
Hard when abundances are super-solar 


e Spectroscopy of Planetary Nebulae: 


Advantages: 


Disadvantages: 


Requires moderate S/N of emission-lines 
Sources easy to find 

Analysis is straightforward 

Can be done in all galaxies 

Measures chemical history of ISM via stars 
Much fainter than H II regions 

Mostly limited to a-elements (best for O) 
Oxygen possibly affected by stellar evolution 
Hard when abundances are super-solar 


e Broadband Photometry of Red Giants: 


Advantages: 


Disadvantages: 


Can be done in all galaxies 

Relatively easy (especially from space) 

Large numbers of stars available 

Contrast is low; difficult to resolve the stars 
Crowding restricts data to low-density regions 
Does not yield data on individual elements 
Mostly limited to estimates of Fe-peak elements 





e X-ray Spectroscopy of Hot Gas: 


Advantages: 
Disadvantages: 


Can be done in distant galaxies/clusters 
Limited (or no) spatial information 
Only measures present day (hot) ISM 
Easy for Fe; difficult for a-elements 
Interpretation can be difficult 


e Absorption-Line Spectroscopy of Background Quasars: 


Advantages: 


Disadvantages: 


Can be done in distant galaxies/clusters 
Can access large array of elements 

Only measures one (halo) sight-line 
Interpretation may be difficult 

Only accesses one ionization state of element 


e Spectroscopy of Integrated Light: 


Advantages: 


Disadvantages: 


Can be done in distant galaxies/clusters 
Measures past ISM via stars 

Can measure both Fe-peak and a-elements 
Requires high S/N aborption-line data 

Only produces luminosity-weighted mean value 
Age-metallicity degeneracy hard to break 
Results are model dependent 








Modeling Chemical Evolution 


To start, let’s consider the types of parameters and variables that 
are involved. First, there are the global variables, all of which are 
a function of time. 


Ez: 


W: 





: Total mass of interstellar gas 


Total mass of stars 


: Total mass of stellar remnants (white dwarfs) 


Total mass of the system 

the rate of mass ejection from stars 
the rate of metal ejection from stars 
the creation rate of stellar remnants. 


Naturally, M; = Ma + Ms + My. 





Next, there are the global parameters of the model, which the 
investigator specifies. Again, all can be a function of time. 


W: Rate of star formation 
f: Rate of infall or outflow of material from the system 
Zr: Metal abundance of the infall (or outflow) material 


(m) : the Initial Mass Function 


As in the equations of photometric evolution, the IMF should be 
normalized so that the total mass is one, as in (8.02). 


Finally, you need four variables which come from stellar evolution 


w: the mass of a stellar remnant 

Tm: the main-sequence lifetime of a star 

Min: the turnoff mass of a population with t = T 

pz: the stellar recyclable mass fraction that is converted to 
metal z and then ejected into space. 





Given the above variables and parameters, the goal is to derive 
Z(t), the fraction of metals (individually, or as a group) in the 
interstellar medium as a function of time. 


Equations of Chemical Evolution 





There are five coupled differential equations which describe the 
chemical evolution of a system. 
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The equations are fairly simple to understand. Equation (9.01) 
is simple mass conservation. Equations (9.02), (9.03), and (9.04) 
keep track of the amount of mass that gets locked up in stars (or 
released into the ISM). Equation (9.05) is the most complex, as it 
describes how the metallicity of the interstellar medium changes 
with time. The first term of (9.05) refers to the amount of ISM 
metals that becomes locked up into stars, the second term gives the 
amount of metals being released by stars, and the third represents 
the amount of metals being brought in (or lost) from outside. 


Although the list of variables mentioned above is formidable, not 
all the variables are independent. Consider E, the mass ejection 
rate from stars. Since mass loss only occurs during post-main 
sequence evolution, the rate of mass loss is related to the number 
of stars evolving off the main sequence at any time. If the galaxy 





consisted only of a single population of stars, this rate would be 
given by 
dMen 

E(t) = Nin(Min — w) = Mo$ (men) (men —w) (8.32) 
However, for a galaxy with on-going star formation, the calculation 
of mass ejection rate must count the main-sequence turnoff stars 
of all stellar ages. Thus the total amount of ISM returned from 
stars at time t is 


pe i Ce ieee EE (9.06) 
where Ma is the upper mass limit of the stellar IMF, and Men, the 
turnoff mass at time t. Similarly, the equation for the total mass 


of remnants formed is 
Mau 


W = w(t — Tm) d(m,t — Tm)dm (9.07) 
Mtn 
The equation for Ez is a bit more complicated since it has two 
terms: one to represent the amount of new metals created by 
a star and released during mass loss, and a second to represent 
the amount of metals that were lost from the ISM when the star 
formed, but are now being re-released. Mathematically, this is 


Ez = Mp-V(t — Tm) (M, t — Tm)dm+ 
[= w= mp.)Z(6 = Fn) = Tn) (Mt = Tm) 


(9.08) 


Finally, there is an equation of metal conservation. If Z, is the 
average metal content in stars, then the total amount of metals 
produced in a galaxy over a Hubble time is 


t Mu 
Z,M,+ZM, = J mp-V(t' —Tm)d(m, t —Tm)dt'dm (9.09) 
0 


Mtn 


Primary and Secondary Elements 


The above equations assume that p,, the fraction of a star which 
is converted into metals, is independent of the initial metallicity 
of the star. In other words, it assumes that the metal under con- 
sideration is a primary element. However, some elements can only 
be made if another element already exists. For example, to make 
nitrogen (via the CNO cycle), the star must already have some 
carbon. Thus, the mass ejection rate of a secondary element, X, 
is 


f= f mpx Z(t E E E t — Tm) dm 


Mtn 


5 (m — w -— mpx)X(t— Tm)¥ (t — Tm) O(m, t — Tm)dm 


Mtn 


(9.10) 


Note that this is similar to the equation for Ez, in that it has two 
terms: the creation term and the recycle term. However, in this 
case, the creation term depends on the prior abundance of Z. 





Analytic Approximation to Chemical Evolution 


Obviously, solving the above coupled differential equations with 
their four free parameters (Y, ¢, f, and Zș) is a non-trivial nu- 
merical problem. However, the problem can be greatly simplified 
if you make two approximations. 


The first approximation to make is to say that the initial mass 
function of stars is independent of time. That is, é(m,t) = o(m). 
Since little is known about how the IMF changes as a function of 
galactic conditions, this may, or may not, be a good assumption. 
(The prevailing theory is that the IMF for a system with no metals 
is heavily biased towards extremely massive stars, but once the 
first bits of metals get introduced into the ISM, this bias goes 
away. ) 


The second approximation is call the instantaneous recycling ap- 
proximation and it is a bit tricker. The approximation says that 
there are two types of stars in a galaxy: those that live forever, and 
those that evolve and die instantaneously. Although this sounds 
like a poor assumption, it’s not as bad as it first appears. Recall 
that the timescales for stellar evolution: 





Main Sequence Lifetimes 


Spectral Mass Luminosity Lifetime 
Type (M/Mo)  (£/Lo) (years) 

O5 V 60 7.9 x 105 5.5 x 105 
BO V 18 5.2 x 104 2.4 x 10° 
B5 V 6 820 5.2 x 107 
AO V 3 54 3.9 x 108 
FO V 1.5 6.5 1.8 x 10° 
GO V 1.1 1.5 5.1 x 10° 
K0 V 0.8 0.42 1.4 x 101° 
MO V 0.5 0.077 4.8 x 101° 
M5 V 0.2 0.011 1.4 x 101 


Note the values. Stars with M > 5Mo evolve in less than 10° 
years, which, in cosmological terms, is almost instantaneously. On 
the other hand, stars with mass less than about 1Mo live forever. 
So the approximation only breaks down for a limited mass range. 


Let’s choose mı to be the dividing line between stars that live 
forever, and stars that evolve instantaneously. Let’s also define 
three new quantities, the Return fraction of gas 


R= ah (m — w)d(m)dm (9.11) 
the Baryonic Dark Matter fraction 
D= wo(m)dm (9.12) 


and the Net Yield (of element i) 


1 
v= TP mp-o(m)dm (9.13) 


M1 
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In words, R is the amount of mass a generation of stars puts back 
into the ISM, D is the amount of mass a generation of stars turns 
into stellar remnants, and y; is the fraction of metal i produced by 
stars for every 1Mo of material locked up into stars or remnants. 
The importance of these three quantities is that each depend only 
on the IMF. If we assume some universal form for the IMF, then 
R, D, and y; are constants that depend only on stellar evolution. 
In other words, they are known quantities. 


Now, let’s take another look at the equations for E, W, and Ez. 
If we assume $(m) is independent of t and use the instantaneous 
recycling approximation, then 


oem f ee E ee en 


Mtn 


= V(t) J “reuco (9.14) 


Mı 


Similarly, the equation for stellar remnants becomes 


wel arten ain 
~ V(t) ‘i ene (9.15) 


and, after a bit of math, 


Ez =W{ZR+y.(1— R)} (9.16) 
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With our two assumptions, the equations of chemical evolution 
become 














a ay (9.17) 

M. = (1— R- D)Y (9.18) 

M = —(1- RY +f (9.19) 

Mu = DV (9.20) 

Ae Me) = —-ZU(1— R)+y.V(1— R) + Zs f (9.21) 


Actually (9.21) can be further simplified by noting that 





d(ZM,) dM 4 dZ 
SS ST, — .22 
dt e Ma dt (922 
Substituting (9.19) for dM4/dt then yields 
dZ 
Mg— = y.¥(1— R) + (Zj — Z)f (9.23) 


For secondary elements, the ejection rate from (9.10) becomes 


Ex = YZ(1 — R)y; + RXV —y.(1- R)XY 
= (1 — R)y-(Z —- X) + YXR (9.24) 
(In many cases, X « Z, so to simplify things, you can often get 


away with Z— X ~ Z.) Metal conservation for secondary elements 
is therefore 


ae) = —XU+U(1— R)(Z — X)ys + RXV + X;f (9.25) 
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But since 


d(XMs)_ydMg,,, dX ga dX 
EE = XE Mg = ~X(1-R)U+X f+ Mg (9.26) 





some of the terms cancel, so 


j= W0- RZ-X) + (X-X) (9.27) 
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Estimating the Stellar Yields 


Models of stellar evolution and supernova nucleosynthesis can pre- 
dict the amount of metals produced by a stellar population. But 
is there a way to empirically check the results of these (purely 
theoretical) calculations? 


Let’s start with the equation for the total amount of metals pro- 
duced over a systems lifetime. Using the instantaneous recycling 
approximation, 


t Mu 
Z3(Ms + Mw) + ZM; = 1 mp. (t' — Tm)o(m, t — Tm)dt'dm 
0 


Mtn 


=| W(t)dt aada 


MY 


=(1-R)y.Vr (9.28) 


where Yr is the total amount of star formation over the history 
of the galaxy. Now that we’ve defined W7, let’s get rid of it. If we 
integrate the equations for the amount of mass that gets locked 
up in metals and remnants, 











t t 
| dMs [4 _R-pD)\w = M,-=(1-R-D)Yr (9.29) 
O 





dt : 
t 
Ee 2G =k Mi we (9.30) 
So 
Z _ (1 — Ry- Yr = M,Z 
5 (1 To R)Vr Ms + Mi 
M,Z 
—— Z —S (eS ee ee . 2 
A= MaE Ma ee 
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Now let’s define u as the gas fraction of a galaxy 


Mg Mg 
KT M+M: + Mo Mi es?) 
So ‘i 
Z, = pa ee -(4)z 9.33 
á Mi — Mu i l—pu l ) 


So, if you choose a region (say, the solar neighborhood), and mea- 
sure the average stellar metalllicity, the present metal abundance 
in the interstellar medium, and the gas fraction, you can estimate 
the stellar yield of the elements in question. 
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Probing Star Formation 


The above equations allow us to estimate the history of star forma- 
tion in the solar neighborhood. For example, if we re-write (9.18) 


as 
1 dM, 


~ T-R_-D dt 


and split the derivative into three parts, we get 


i= 1 (e) (25) ( dM, ) (9.35) 
1-R-D\ dlogZ dt dlog Ms 

Now consider the derivatives. Since we can measure the metallicity 
of solar neighborhood stars, we can determining how much stellar 
mass there is as a function of metallicity. Thus the first derivative 
in (9.35) is a measureable quantity. Similarly, if we study nearby 
F-stars, and compare their absolute luminosities to the F-star zero- 
age main sequence luminosity, we can estimate how much main- 
sequence evolution has occurred, i.e., we can estimate their ages. If 
we measure the stars’ metallicities as well, then we have dlog Z/dt. 
Finally, the last derivative is simply In(10)M,. Thus, we can 
measure the history of star formation in the solar neighborhood. 





wv (9.34) 





This leads to the long-standing G-dwarf problem. Simple models 
of chemical evolution (say, a closed-box model or a constant in- 
fall model) predict many more low-metallicity stars in the solar 
neighborhood than are observed. 
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Figure 4. Left: Local Age-Metallicity relationship, in 1 Gyr age bins (top) and in 0.1 dex metal- 
licity bins (bottom); isocontours correspond to the data of GKS (Nordstrom et al. 2004), large 
symbols show the corresponding relationships of age vs mean metallicity (top) and Metallicity 
vs mean age (bottom), small symbols the complementary relationship, and data with error bars 
designate results of independant studies with different samples (from Subiran et al. 2007 and 
Subiran and Girard 2005). Notice that in both panels the Sun appears to be slightly metal-rich 
for its age or slightly older than stars of the same metallicity. Right: By using simulated data 
as input for the AMR (points with 1-o Gaussian error bars in lower panel, as indicated by dotted 
lines) and a star formation rate exponentially increasing with time (top panel), one may derive 
the corresponding age vs average metallicity (circles) and metallicity vs average age (squares) 
relations; the former is flatter than and the latter is steeper than the input relationship (as with 
the case of the real data, see left pannels). Their average (solid curve) is a better approximation 
to the input data. 
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[Casuso et al. 





To get the history of matter infall, we can do a similar manipula- 
tion with equation (9.19) 
dM 
= 1— R)W 
a ) 
_ dM, 1-R dM, ) dlogZ 
= dlogZ 1—R-—DdlogZ dt 














(9.36) 


If we measure the metallicities of clouds of H I and Hə of different 
masses, then all terms of this equation are known, and the history 
of matter infall can be found. 
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The Closed Box Model of Chemical Evolution 


As an example of what a chemical evolution model can do, consider 
a closed system, where all the material for current star formation 
comes from mass lost by a previous generation of stars. In this 
case, there is no infall, and, from (9.23), 


dZ 


Mg— = yeV(1—R)+(Zp—Z)f=yeV(l-R) (9.37) 


In addition, from (9.19), we have 


we ae Ge 2) eae eee GN ca (9.38) 





By dividing these two equations, we get 


dZ /dM, dZ 
—= — — >=> — Z 9.39 
Mend a aM ee) 


Since y; is a constant of stellar evolution 


Z M 
d 1 dM M 
dZ = -y. | z = Zi - o= —Yz In (=) 


(9.40) 
where Z and M,, represent the initial metallicity and gas mass 
of the galaxy, and Z; and M,, represent those quantities today. 


Now let 
= —_—_ = — 41 
j (Z) 7 (2) ô (22) een 


Time ln (=) (9.42) 
Ho 











SO 
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or 
Hı = Ho exp j- a-a) (9.43) 
Yz 
In other words, as the system evolves, the gas fraction will decrease 
exponentially with Z. Of course, we can’t measure the metallicity 
evolution of the ISM directly, but we can use stars as a probe. If 
we take the derivative of (9.43) with respect to Z, then 


du H1 Zı — Zo 
eee eae 44 
dZ Yz mer Yz eae 








Meanwhile, through(9.18) and (9.19) 


dM, [/dM, do (1—R-D) 





| Ge= Gs (9.45) 
(2) (5) -(2) 2) ol 54) 
(9.46) 


Finally, if we put this equation in terms of log Z, instead of Z, 
then 


le mo (2) (BGP) (2) {222 
dlog Z Us 1—R O71 Ye 

(9.47) 
Because the number of stellar metallicity measurements has (tra- 
ditionally) not been extremely large, many times people plot the 
cumulative distribution, i.e., the number of stars with metallici- 
ties less than Z. This is simply found by integrating (9.47). If we 
collect the constant terms and let 


(A) e 
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then for the cumulative distribution 


aS L-G fexp (-4-*) -1} (9.49) 
01 Yz 


When this is fit against observations, it is clear that the metallic- 
ities of stars in the solar neighborhood cannot be fit with a closed 
box model of galactic evolution. Either f Æ 0, or Zo Æ 0, or 
the initial population of stars did not have the same IMF as the 
stars today, or there are severe chemical inhomogeneities in the 
ISM, and star formation occurs preferentially in regions with high 
metallicity. 





A similar closed-box calculation can be performed for secondary 
elements. For these, if you divide (9.27) by (9.19), then 








dX /[dM, dX Wi1-— R)\(Z— X)ye 
ee, = ee ee 9.50 
me dt gee dM, —WU(1— R) (20) 
If we assume that X < Z, then this simply reduces to 
M,dX = —y,ZdM, (9.51) 
Now if we assume that M, = M; at t = 0, then from (9.42) 
Z,—-Z 
Myg = Mi exp {= oI (9.52) 


or 





m=- (4) cp { 2-2} az = (Me) az os 


Thus i 
M,dX = “7247 => X= 5 (=) z? (9.54) 
Yz z 
In other words, if X is a secondary element, then a plot of log- 
X versus log-Z (i.e., [X] vs. [Z]) should have a slope of 2 (and 


presumably go through the solar value). 
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Timescales for Chemical Evolution 


The timescales for chemical evolution are simple to derive. First, 
let’s consider the gas consumption timescale. In the absense of 
accretion, how long does it take to use up the gas? 


MyM, 
nam, / Hel- apr a R O 


In the solar neighborhood, M, ~ 5.7Mo pc™?, and the star for- 
mation rate is VU ~ 4.2Mọo pe? Gyr”. For R ~ 0.4, that means 
that the gas should e-fold in Tą ~ 2 Gyr, much less than a Hubble 
time. 


Similarly, we can calculate the timescale for chemical enrichment. 


zl \ el 
M Z 


4 Tal 
= ee 9.56 
ER a 0) 


Since y ~ Z, Tz ~ Tą. In other words, for the solar neighbor- 
hood, the metallicity (and amount) of the gas should e-fold rather 
quickly. This strongly suggests that some our region of space has 
received matter from some other region. 
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The Balanced Infall Model 


A closed-box model is not realistic for most systems: there is good 
evidence that infall of new material plays a part in the chemical 
evolution of a region. So let’s calculate (for simplicity) a model 
where the infall just balances the amount of material becoming 
locked up into stars. In other words, a model where the mass in 
the interstellar medium stays constant. In this case 


dMy 
dt 





=-V+E+f=—-(1-R)V+f=0 (9.57) 


So 


oe i 
v= (9.58) 


The equation for metals is then 


M, =W- Ry: + (Zp —Z)f = flue + Z-Z) (9:59) 


Also, since there is accretion 











dM, _ 
f (9.60) 
So 
dZ | dM, fete DZ 
— = Geen Ooae Z A a Z 
Tat / dt IdM; f errar 
(9.61) 


If we integrate this 


Z M 
1 dZ d 
1 —— = a (9.62) 
Zo Yet Ze-Z Mo Mg 
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then 


Inj BEE eee (9.63) 


Ue Le HLA Mg 


Now let v represent the total amount of mass accreted, scaled to 
the mass in the ISM, i.e., v = (M —Mo)/M,. Then 


Zi = (y + Zp — e`”) + Zoe” (9.64) 
If the galaxy began with Zo ~ 0 and Zş ~ 0, then 
Zi = y-(1 = e) (9.65) 


Note that if you also assume that the galaxy began as an entirely 
gaseous system, then 


M-Mo M-My 3 
— Se OSOS —_| 1 . 
E M; Mo d (200) 


Thus, as u — 0, v — œ, and 4, — y,. In other words, the 
metallicity of the system asymptotes out at the value of the stellar 
yield. 


The balanced infall model is similar to the closed-box model in 
that Z x yz, so the metallicity of a system is actually measuring 
the stellar yields. Also, note that Z/y, is not a strong function 
of the gas fraction, or V/f. It depends mostly on the current 
properties of the system. 
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The Mass-Metallicity Relation for Galaxies 


The chemical evolution formulism described above provides a con- 
text to explain the correlation between luminosity (mass) and 
metallicity seen in galaxies. For simplicity, let’s consider an el- 
liptical galaxy. As you recall, there is a correlation between the 
absolute magnitude of ellipticals and their color, with large bright 
galaxies being redder. This is often interpreted as a metallicity 
effect. 


From (9.09), the total mass of metals formed through the history 
of a galaxy is 


t OO 
ZM: = J Í mp-V(t — Tm)d(m, t — Tm)dt'dm (9.67) 
0 


my! 


which, in the instantaneous recycling approximation, is 


t oO 
ZM: =j W(t)dt mpzo(m)dm (9.68) 
O My 
From (9.28) and (9.13), this is 
Mes 1- R 


ZM: Gad R) =y.M (9.69) 


hee {TRE 


In the early stages of galaxy formation Ma ~ M;, and D ~ 0, so 


yzM, 
ARN A ; 
M (ARU) 


t 


Now we can use two possible arguments. First, let’s suppose that 
the supernova rate early in the galaxy’s history was high, and 
that these explosions pumped enough energy into the interstellar 


24 


medium to cause the gas to escape. Since the number of super- 
novae is proportional to the number of stars in the galaxy 


GMM, GM? M? 


SE 2L .71 
R R sN X M, => R x Ms (9.71) 





Now, if we assume some mass-radius relation, M+ x R°, then 
Mi 

R 
But from (9.70), M; x (M,/Z), so 


x MI) x Ms (9.72) 


M 2—(1/a) ea 
(=) x Ms => Z x M2 (9.73) 


To get the (very rough) empirical value of Z ~ M?:?5, œa ~ 1.5. 


Alternatively, suppose that elliptical galaxies were formed during 
proto-galaxy collisions. If the efficiency of star formation depends 
on the mass of the system 





Ww Ms 
(5) x MP —> (2) x MiM? = M, x MPTE 
(9.74) 


In the above proportionalities, we’ve assumed that early on, Mg ~ 
M+, and the mass of stars we see today, M,, is proportional to 
the original star formation rate. With this relation, and (9.70), we 
get 

M 


p+1 p 
Ze) — Zx Mg” (9.75) 


p+1 
Ms x M4 x ( Z 


and p ~ 1/3 recovers the empirical relation. Note that in both 
cases, Zs — yz, So eventually the mass-metallicity relation has to 


flatten out. 
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